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Abstract: We derive universal scaling laws for physical observables such as the critical 
temperature, the chiral condensate, and the pion decay constant as a function of the 
flavor number near the conformal window of many-flavor QCD in the chiral limit. We 
argue on general grounds that the associated critical exponents are all interrelated and can 
be determined from the critical exponent of the running gauge coupling at the Caswell- 
Banks-Zaks infrared fixed point. We illustrate our findings with the aid of nonperturbative 
functional Renormalization Group (RG) calculations and low-energy QCD models. 
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1. Introduction 

Many-flavor QCD has recently attracted a great deal of attention for a number of reasons: 
first, as a controllable deformation of real QCD, it can teach important lessons about the 
chiral structure of QCD-like theories. Second, it serves as a building block for alternative 
technicolor- like scenarios for the Higgs sector. And third, it exhibits a quantum phase 
transition from the QCD-like chirally broken to a conformal phase as a function of the 
flavor number and thus gives rise to interesting quantum critical behavior. 

Adding A'^f massless quark flavors to a nonabelian SU(A^c) gauge theory increases 
the screening property of fermionic fluctuations. An obvious consequence is the loss of 
asymptotic freedom for Af > A^f^'^' := ^Ac (= 16.5 for SU(3)). Already at smaller Af, 
A'f > i^fJi'L^ (— 8.05 for SU(3)), the second (3 function coefficient changes sign induc- 
ing an infrared (IR) attractive fixed point of the gauge coupling a* > [Q. For small 
jya f- — _/Vf > 0, the fixed-point value is small. This gives access to a perturbative analysis, 
suggesting that the system approaches a conformally invariant limit in the deep IR . For 
smaller A^f (such as real QCD), this Caswell-Banks-Zaks fixed point is destabilized due to 
the spontaneous break-down of chiral symmetry, resulting in massive fermionic excitations, 
strongly-coupled glue and massless Goldstone bosons. 

The above considerations propose the existence of a critical fiavor number A^^'^'^, sep- 
arating the chiral-symmetry-broken phase from the conformal phase. Theories with a 
flavor number A^f satisfying A^^'^'' < Af < A^^^'^' are said to be in the conformal window. 
The fiavor number A^f therefore serves as a control parameter for a quantum phase tran- 
sition. Investigations of this phase structure have been performed by continuum meth- 
ods H, I, I, I, 0, I, |, In], H, m, m, 0], as well as lattice simulations 



n, m, m, m, m, H, m, H, l^, H, H, H, Hl- Recent results have collected a 
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substantial body of evidence for the existence of the conformal phase and a critical fla- 
vor number, defining the onset of the conformal phase: lattice simulations have provided 
evidence that 8 < A^"" < 12 |5|, 1^, ||, ||], even though the case Nf = 12 is contro- 
versial, see |31]. On the present level of accuracy, these results go well together with an 
earlier quantitative estimate from a combination of four-loop perturbation theory with the 
functional Renormalization Group (RG) which yields A'^^'^'' = 10.0 ['[i'^ [12|. 



Beyond the precise location of this quantum critical point on the A^f axis, the critical 
behavior in the vicinity of the fixed point is expected to show several peculiarities: for 
A'^f < A'^f'^'' in the chirally broken phase, a standard chiral phase transition can be antic- 
ipated to occur at finite temperature. The critical phenomena near this finite-T chiral 
phase transition are expected to be determined by the precise symmetries of the order 
parameter, i.e., the chiral condensate, defining the universality class of the transition (e.g., 
SU(2)l(8) SU(2)r ~0(4) for A^f = 2) |Q. Near the critical temperature an effective descrip- 
tion in terms of a Ginzburg-Landau-type effective potential is expected to hold even on a 
quantitative level. By contrast, the zero-temperature quantum phase transition as a func- 
tion of the control parameter A'f may not have a continuous Ginzburg-Landau description. 
In particular, there appear to be no light scalar states in terms of which an effective theory 
could be constructed on the conformal side A'^f > A^^" of the phase transition [^, 
even though the order parameter (chiral condensate) should change continuously across 
the phase transition. 



Plotting the chiral-phase-transition temperature versus the control parameter A^f yields 
the phase boundary in the (T, A'f) plane which has first been explored in [13, |l^ using the 
functional RG. Near the critical flavor number, an intriguing relation between the shape 
of the phase boundary and RG properties of the Caswell-Banks-Zaks fixed point has been 
identified. This relation connects the scaling of the critical temperature with the critical 
exponent of the gauge coupling near the Caswell-Banks-Zaks fixed point. In the present 
work, we argue that this scaling relation can be extended to further physical observables 
such as the chiral condensate or the pion decay constant. As these observables are accessible 
to a variety of nonperturbative methods, our results suggest that the corresponding scaling 
can become a useful tool to study the phase transition to the conformal phase quantitatively. 



Our main arguments are based on very general considerations and involve only few 
assumptions about the RG structure of the theory. These are illustrated in Sect. ^ for the 
simple few-flavor case and in Sect. ^ for many flavors near the conformal window. These 
arguments are then made more concrete with the aid of functional RG calculations in a 
derivative expansion of the QCD effective action, or simply within low-energy QCD models 
in Sect. ^. Throughout this work, we concentrate on the chiral phase transition even though 
we also expect an impact of the confining nature of the theory on the properties of the 
system near criticality. However, as we work in the chiral limit, there is no good order 
parameter for confinement, implying that nonanalyticities in the correlations are rather 
dominated by the chiral degrees of freedom. 
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2. Few-flavor QCD and the role of scale fixing 

QCD in the chiral limit of zero quark masses depends only on one parameter in the Eu- 
clidean Lagrangian, namely the gauge coupling 5, 

^ = i^^M-^;- + i^^t^^MM^- (2-1) 

In the quantum theory, the gauge coupling has to be fixed at a certain momentum scale in 
terms of a renormalization condition. The renormalization group finally trades the gauge 
coupling fixed at an arbitrary scale for one single parameter Aqcd of mass dimension one 
which sets the mass scale for all physical obscrvablcs of the theory. In other words, all 
physical observables respond trivially to a variation of Aqcd according to their engineering 
dimension. In units of Aqcdj the theory is completely fixed. 

In order to discuss the dependence on quantities such as the fiavor number, it is 
important to emphasize that a variation of the flavor number docs not merely correspond 
to a change of a parameter of the theory. It rather corresponds to changing the theory 
itself. In particular, there is no unique way to unambiguously compare theories of different 
flavor number with each other, as different theories may have different scales Aqcd- 

For instance, it might seem natural to compare theories with different flavor numbers 
at fixed Aqcd with each other. But Aqcd itself is not a direct observable, so that such 
a comparison is generically inflicted with theoretical uncertainties. Moreover, Aqcd is 
regularization-scheme dependent which can affect comparisons between different theoret- 
ical methods, say, lattice and continuum results. Another option could be a scale flxing 
in the deep perturbative region, say, at the Z mass pole by fixing a{Mz)- However, theo- 
ries with different flavor numbers then exhibit a different perturbative running, such that 
IR observables vary because of both high-scale perturbative as well as non-perturbative 
evolution. 

Instead, wc propose to choose a mid-momentum scale for the scale fixing, as the high- 
scale perturbative running is then separated from the more interesting non-perturbative 
dynamics. In this work, we fix the theories at any A^f by keeping the running coupling at the 
T mass scale fixed to a{mr) = 0.322. Even though also this choice is scheme dependent, 
these dependencies should be subdominant, as they follow a perturbative ordering. In 
general, fixing the scale via the coupling is a prescription which is well accessible by many 
nonperturbative methods. 

Let us now present a simple argument that illustrates how A^f dependencies of physical 
observables can be understood in the limit of small Nf. As already stated above, all IR 
observables such as the critical temperature Tcr, pion decay constant f-j^, chiral condensate 
(■j/i'i/;)^/^, and model-dependent concepts such as the constituent quark mass, are propor- 
tional to Aqcd- The latter on the one hand can be read off from the UV behavior of 
the running coupling, a{k) l/ln(fc/AQCD) for large k. On the other hand, the value of 
Aqcd can be associated with the position of the Landau pole in perturbation theory.^ In 

^Of course, this statement has to be taken with care, since Aqcd is a meaningful scale, whereas the 
Landau pole is simply an artifact of perturbation theory. 
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this simple sense, the artificial Landau pole can be taken as an estimate for the scaling 
of physical observables. In one-loop RG-improved perturbation theory, the position of the 
Landau pole can be read off from 

0^—^='+4.bo In^, 
a(AQCD) a{^o) /^o 

^ ^^A^.-?iVfl, (2.2) 



where fiQ denotes a perturbative scale, such as m-r, Mz, .... Solving this equation for Aqcd 
and expanding the result for small A'^f leads us to 



Aqcd — Moe ^Trboo^c^o) 

~ Ho e" iiiVc^(Mo) (1 - eiVf + ©((eiVf)^)) . 

Choosing fiQ = ""^r, we find e = x2ijv|a(^o) ^ 0.107 for Nc = 3. Two conclusions can imme- 
diately be drawn: first, Aqcd can be expanded in A'^f and has a generically nonvanishing 
linear term; second, for the present way of scale fixing, the linear behavior should be a 
reasonable approximation for finite values of A^f, say A'f < 4, as the expansion parameter e 
is small. 

As Aqcd sets the scale for all dimensionful IR observables, we are tempted to conclude 
that all IR observables scale linearly with A'^f for small A'f with the same proportionality 
constant e. This is, of course, a bit too simplistic, as the dynamics which establishes the 
value of the IR observables generically carries an A'f dependence as well. E.g., the chiral 
symmetry-breaking dynamics depends on the number of light mesonic degrees of freedom, 
which is an A'f-dependent quantity. Detailed quantitative model studies ||3^, however, 
demonstrate that the scaling of the critical temperature does not receive strong corrections 
of this type and indeed scales according to 

T,, = To{l-eNi + 0{{eNif)), (2.3) 

where Tq is a dimensionful proportionality constant. We conclude that the phase boundary 
in the (T, A'f) plane has a linear shape for small A'^f which can mainly be understood as a 
result of the perturbative A'f scaling of Aqcd- Note that this observation is consistent with 
lattice simulations |36] and has been exploited for parameter fixing in PNJL/PQM-model 



studies 1 37 ] . In the following, we will show that the shape of the phase boundary as well as 
the A'^f scaling of other observables for large A^f can also be understood from simple scaling 
arguments which this time follow from general properties of the nonperturbative domain. 

3. Scaling in many-flavor QCD near the conformal window 



In this section, we review, detail and extend the scaling arguments presented in |13, 14 1, 
leading to universal relations near the conformal window. Whereas the A'^f scaling in the 
few-fiavor case essentially follows from analyticity of the observables in A^f also for A'f near 
zero, the A'f dependence near the conformal window is clearly nonperturbative in A^f. The 
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Figure 1: Left panel: illustration of the IR running of the gauge coupling in comparison with the 
critical coupling g^^ for xSB. Below the conformal window, iVf < N" , exceeds the critical value, 
triggering the approach to xSB. For increasing flavor number, the IR fixed-point value drops 
below the critical value denoting the onset of the conformal window. Right panel: sketch of the /3 
function of the running coupling. The slope of the P function at the IR fixed point corresponds to 
minus the critical exponent 0, cf. Eq. (3.1). 



mere existence of the conformal window requires a sufficient amount of fermionic screening, 
i.e., a sufficient amount of fermionic degrees of freedom A'^f > A^^'^''. 

The lower end of the conformal window is characterized by the onset of chiral symmetry 
breaking. Whereas the coupling approaches the Casweh-Banks-Zaks fixed point gl in the 
conformal window, chiral symmetry breaking destabilizes this fixed point below the critical 
flavor number A'^^'^''. This suggests the existence of a critical coupling g1^. If > g^^ at 
some scale, the system will be triggered to run into the chirally broken regime. 

For a monotonous coupling flow, the value of the Caswell-Banks-Zaks fixed point gl or 
its nonperturbative variant corresponds to the maximum possible coupling strength of the 
system in the conformal window.^ As both gl and g^^. depend on the number of flavors, the 
condition g1{N^^) = g'^j.{N^^) defines the lower end of the conformal window and thus the 
critical flavor number, see left panel of Fig. ^ for an illustration. For g"^ > g"^^, the model is 
below the conformal window and runs into the broken phase. Slightly below the conformal 
window, the running coupling g^ exceeds the critical value while it is in the attractive 
domain of the IR fixed point g"^. The flow in this flxed-point regime can approximately be 
described by the linearized /3 function 

/3g2 ^ dtg^ = -e {g^ - gl) + - glf), (3.1) 

^The following scenario does not apply if the running coupling overshoots, develops a local maximum and 
approaches the fixed point from above which requires a double- valued [3 function. This behavior is seen, e.g., 
in MOM-scheme running couplings derived from the ghost-gluon vertex using truncated Dyson-Schwinger 
equations for pure Yang-Mills theory. If such a scenario held in many-flavor QCD, the occurrence of 
XSB, i.e., whether or not the fermion sector becomes critical for a given Nf, would depend quantitatively on 
the details of the coupling flow, making it difficult to extract universal features. However, we do not expect 
such a running to occur in many-flavor QCD near the conformal window as the fast running of fermions 
in the near-critical region reduces the fermionic screening contributions, thus supporting a monotonous 
increase of the gauge coupling. Also, such a behavior is not observed near the upper end of the conformal 
window where perturbation theory is expected to hold. 
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where t = ln(A;/A) with k being a suitable RG scale and A defining the UV cutoff. The 
universal "critical exponent" Q denotes (minus) the first expansion coefficient. We know 
that < 0, since the fixed point is IR attractive, see right panel of Fig. |l[ In general, the 
critical exponent depends on A^f, = 0(A^f). The solution to Eq. ( |3.1] ) for the running 
coupling in the fixed-point regime reads 

where the scale ko is implicitly defined by a suitable initial condition and is kept fixed in the 
following. It corresponds to a scale where the system is already in the fixed-point regime, 
and otherwise plays the same role as the renormalization scale /Uq in Eq. ( |2.2D ; in particular, 
physical observables are independent of ko- For the present fixed-point considerations, it 
provides for all dimensionful scales in the following. But knowing the full RG trajectory, 
ko is related to fio and thus, say, to the initial r mass scale by RG evolution. 

Our criterion for xSB is that g'^ik) should exceed for some value of A; < kcr- From 
Eq. ( |3.2[ ) and the condition ^^(/ccr) = '^^ derive the estimate valid in the fixed-point 
regime 

kcr^ko{gl-gi)~^. (3.3) 

This scale fccr now takes over the role of the fixed renormalization scale = m,- in the 
small- A^f argument given above: it sets the scale for the critical temperature Tcr ~ /ccr with 
a proportionality coefficient provided by the solution of the full flow. The last step of the 
argument goes along with the estimate that the IR fixed-point value gl roughly depends 
linearly on A'^f. More precisely, we assume that the A'f dependence of the coupling quantities 



can be linearized near the critical flavor number. From Eq. (3^), we thus flnd the relation 



T,,^ko\Nf-N^'\-B, (3.4) 

which is expected to hold near A'^f'^'' for A'^f < N^^. Here, & should be evaluated at N^^."^ 
Relation ( p.4| ) is an analytic prediction for the shape of the chiral phase boundary in the 
{T,N{) plane of QCD. This result is remarkable for a number of reasons: first, it relates 
two universal quantities with each other: the phase boundary and the IR critical exponent. 
Second, it establishes a quantitative connection between the chiral structure (Tcr) and the 
IR gauge dynamics (0). Third, it is a parameter-free prediction following essentially from 
scaling arguments. 

As Eq. (|3.4| ) relates two universal quantities, it is important to understand to what 
extent the underlying argument makes use of non-universal but scheme-dependent quanti- 
ties. Of course, the running coupling is a strongly scheme-dependent concept, and so is the 
value of the IR fixed point g"^. However, the existence of the fixed point as well as the value 
of the critical exponent are scheme independent.^ Also, the value of the critical coupling 



'Accounting for the A'^f dependence of O by an expansion around A'f^'' yields mild logarithmic corrections 



to Eq. ( |3.4 - 

*The running coupling in addition is definition dependent. Here, we assume that the running coupling 
used in the discussion provides for a reasonable measure of the interaction between the gauge and the quark 
sector which manifestly exhibits the IR fixed point in the conformal window. 
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fl'^j, is scheme dependent. Nevertheless this scheme dependence has to cancel against that 
of the running coupling itself, as whether or not xSB occurs is a universal feature of the 
system. 

Furthermore, we have implicitly neglected all dependencies of the running on scales 
other than the RG scale k. At finite temperature, there will, of course, be dependencies of 
the coupling on T in particular at low scales where T/k becomes large. The relevant scale 
for the above argument, however, is the scale kcr where xSB is triggered. We expect that 
this scale is generically somewhat larger than the temperature as long as we are in the xSB 
phase or approach the phase transition from below. Indeed, this assumption turns out to 
hold in all model calculations as well as in the RG results described below. This allows us 
to ignore the T dependence of the running coupling g'^ and of the critical coupling gcr- 

Let us now generalize these considerations to other physical observables such as the 
chiral condensate or the pion decay constant. To be more specific, we are interested in 
the scaling behavior of physical observables as a function of the number of massless quark 
flavors A^f at vanishing temperature. In this case, the above argument can be followed 
straightforwardly, where the determination of the critical (RG) scale kcr in Eq. ( |3.3D plays 
a prominent role. In terms of low-energy effective theories, the scale kcr can be viewed 



as an ultraviolet (UV) cutoff. From Eqs. (3^) and (3^), it follows immediately that the 



critical scale and therewith this UV cutoff Agfj of the low-energy sector is tightly related 
to the critical flavor number: 



Aes^kcr^ko\N{-N^'\-B. (3.5) 

On the other hand, observables O with mass dimension dcj which are computable in this 
effective field theory defined by the fixed-point regime are necessarily related to the UV 
cutoff Aeff in a simple manner, 

O coAt§ , (3.6) 

where co is a numerical constant which depends on the details of the theory, e.g., the 
number of colors A^c and also the number of flavors Af. Here, we have assumed that a 
general separation of scales holds in the sense that all effective UV parameters of the low- 
energy effective theory are fully determined by the quark-gluon dynamics in the mid- and 
high-momentum regime. Combining Eqs. (3.5) and ( |3.6| ), we find 



O ^k^^'lNi- N^'\—8-. (3.7) 

This relation extends the scaling properties of the critical temperature near the conformal 
window found above in Eq. (3.4) to that of other physical IR observables. Again the 



universal scaling of these observables as a function of Af is related to the critical exponent 
& of the running coupling. 

4. Functional RG results 

The critical temperature as a function of A^f has first been computed in 14] in the 



framework of the functional RG |3^ (see |40] for reviews on the functional RG in gauge 
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theories). In this section, we briefly review these results and discuss them in the Ught of 
the scahng relations. Result on further IR observables will be discussed in the next section. 

In H, H], the RG flow of QCD starting from the microscopic degrees of freedom 
in terms of quarks and gluons was studied within a covariant derivative expansion, ap- 
proaching the critical xSB temperature from above. A crucial ingredient for xSB are the 
scale-dependent gluon-induced quark self-interactions of the type 



rv>,mt = J Kp'yS'^ai^ptp-ytpS, (4.1) 

where a, . . . denote collective indices including color, flavor, and Dirac structures. These 
four-fermion interactions are set to zero at the initial UV scale, \ai3'ys\k^A ~^ 0. This 
guarantees that the A's at A: < A are solely generated by quark-gluon dynamics from first 
principles (e.g., by IPI "box" diagrams with 2-gluon exchange). This is an important 
difference to models such as the Nambu-Jona-Lasinio (NJL) model, where the A's are 
independent input parameters. 

As we approach the chiral phase transition temperature from above, the derivative 
expansion with local "point-like" interactions is a self-consistent approximation. This cor- 
responds to replacing the momentum structure of the A's by the overall k dependence 



on the RG scale. This approximation has been successfully tested quantitatively in |12| 
by verifying the resulting insensitivity of the many-flavor quantum phase transition on 
the momentum regularization. In the chirally broken regime, this approximation breaks 
down as, e.g., mesons manifest themselves as momentum singularities in these vertices. 
These restrictions result in a total number of four linearly independent couplings Aj |41|. 



Introducing the dimensionless couplings Aj = /c^Aj, the corresponding /3 functions read 

dtXi = 2Xi - bijXjO^ - AijkXjXk - Cig'^, (4.2) 

where the coefficients A, b, c depend on the temperature, number of quark flavors Nf and 
number of colors N^, for explicit representations, see |l^, p^ . 

Within this truncation, a simple picture for the chiral dynamics arises, see Fig. |2|: 
at weak gauge coupling, the RG flow generates quark self-interactions of order A ~ 
via the last term in Eq. (|4.2| ) with a negligible back-reaction on the gluonic RG flow. If 
the gauge coupling in the IR remains smaller than a critical value g < gcr, the A self- 
interactions remain bounded, approaching flxed points A* in the IR. The fixed points 
correspond to a shifted Gaufiian flxed point A^^"'^ 1^2=0 = 0. At these flxed points, the 
fermionic subsystem remains in the chirally invariant phase which is indeed realized at 
high temperatures T > Tcr- 

The evaluation of the QCD RG flow in a covariant derivative expansion |42, 43 1 in- 



cludes further gauge-fleld operators as well as kinetic terms for the fermion. From the gauge 
dynamics, the running gauge coupling can be extracted including its dependence on tem- 



perature T as well as flavor and color numbers. Following Ref. |13, 14], the increase of the 
running coupling in the IR is weakened on average for both larger T and larger A^f, in agree- 
ment with general expectations. In addition, also gcr depends on T and Nf, even though the 
A^f dependence is rather weak. For instance, the (non-universal) zero-temperature value of 
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Figure 2: Sketch of a typical /3 function for the fermionic self-interactions (taken from [p^): 
at zero gauge coupling, g — (upper solid curve), the Gaufiian fixed point Ai = is IR attractive. 
For small g ^ (middle/blue solid curve), the fixed-point positions are shifted by the gauge- field 
fluctuations ~ g'^. For gauge couplings larger than the critical coupling g > g„ (lower/green 
solid curve), no fixed points remain and the self-interactions quickly grow large, signaling x^B. 
For increasing temperature, the parabolas become broader and higher, owing to thermal fermion 
masses; this is indicated by the dashed/red line. 



the critical coupling for an optimized RG scheme is acr = glr/i^'^) ~ 0-^ A^c = 3 and a 
wide range of A^'f jl^ . 

The T dependence of g^r arises from the quark modes acquiring thermal masses. This 
leads to a quark decoupling, requiring stronger interactions for critical quark dynamics. In 
the /3 function picture of Fig. ^, the Aj parabolas become broader with a higher maximum; 
hence, the annihilation of the Gaufiian fixed point by pushing the parabola below the Aj 
axis requires a larger gcr- 

At zero temperature and for small A'f, the IR fixed point g1 is far larger than g'^^, 
hence QCD is in the xSB phase. For increasing T, the temperature dependence of the 
coupling and that of g^-^ compete with each other. For the case of many massless quark 
flavors A'f, the critical temperature is plotted in Fig. |3[ For the scale fixing at the r mass 
scale discussed above, we observe an almost linear decrease of the critical temperature for 
increasing A'^f with a slope of AT^ = T{Nf) — T{Nf+ 1) ~ 25MeV at small A'^f. This linear 
dependence of the full result confirms the simple estimate given in Eq. ( p.3| ) to a very good 
accuracy. Moreover, the predicted relative difference for Tcr for A'^f = 2 and 3 flavors of 
A ~ 0.146 is in very good agreement also with lattice studies [36|. We conclude that the 
shape of the phase boundary for small A'^f is basically dominated by fermionic screening. 

For larger flavor number, the critical temperature decreases and the phase transition 
line terminates at the zero-temperature quantum phase transition at N^^, denoting the 



onset of the conformal window. In the RG study of |T^, 14], we flnd a critical number of 
quark flavors, N^'^ ~ 12.9. This result for N^'^ agrees with other studies based on the 2-loop 
/3 function However, the precise value of A'^^'^'' has to be taken with care: for instance, 
in a perturbative framework, A^^'^'" turns out to be sensitive to the 3-loop coefficient which 
is not reliably reproduced in this leading-order study. This coefficient can bring A'^^'^'' down 
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Figure 3: Chiral-phase-transition temperature T^r versus the number of massless quark flavors Nf 
for N[ > 2, as obtained in Ref. The flattening at A^f > 10 is a consequence of the IR fixed-point 
structure. The dotted line depicts the analytic estimate near iV"' which follows from the fixed-point 
scenario (cf. Eq. (3.4)). 



10.0+0^ which remains stable under the inclusion of 4-loop corrections |T^. The 



to 

error bars parameterize truncation errors which habe been quantified by artificial scheme 
dependencies. 

To the left of the conformal window (A^f < N^^), the phase transition line shows a 
characteristic flattening. This is again in perfect agreement with our scaling relation (^.4|). 
The fit to numerical results from a functional RG approach is depicted by the dotted line 
in Fig. ^. In particular, the fact that |0| < 1 near N^^ explains the flattening of the phase 
boundary near the critical flavor number. Within the covariant derivative expansion of RG 
flow, the IR critical exponent of the /3 function of the coupling at the critical flavor number 
yields Q{N^'^) ~ —0.60 [14|. However, this estimate is likely to be affected by truncation 
errors as Q{Nf) is expected to show sizable dependencies on A'^f; hence, any error in N^'^ 
translates into a corresponding error in Q{N^^). 

For comparison, we plot the perturbative estimates for the critical exponent Q{N{) in 
Fig. ^ based on the 4-loop /? function in the MS scheme [^] . The apparent convergence of 
the perturbative expansion is remarkable as the difference between the 3- and 4-loop result 
is below the 1% level (for 8.5 < A'f < 16.5). The 2-loop result shows larger deviations for 
smaller N{, as the fixed-point coupling is larger in this regime. In Fig. ^, we also show 
the estimate of from the covariant derivative expansion of the RG flow |14]. While this 
estimate includes nonperturbative contributions to all-loop orders in the gauge sector, the 
derivative expansion in the fermion sector effectively corresponds to the inclusion of just 
the (RG-improved) one-loop quark diagram. This explains a large part of the difference to 
the perturbative estimates. Incidentally, results from further studies of the full /3 function 
can be used to estimate 0. For instance, the values obtained from the conjectured 



"NSVZ-inspired" /3 function [45] are identical to the 2-loop result. 
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10 12 14 16 18 

Nf 

Figure 4: Critical exponent \Q\ of the running gauge coupling at the Caswell-Banks-Zaks fixed 
point as a function of the flavor number N{. The perturbative expansion appears to converge, as 
the difference between the 3-loop (blue/dashed line) and 4-loop result (red/solid line) is below the 
1% level (for 8.5 < Nf < 16.5). The 2-loop result (green/dot-dashed line) shows larger deviations 
for smaller Nf. The black/solid line with symbols corresponds to the estimate from the covariant 
derivative expansion of the RG flow . The inlays depict the characteristic shapes of the scaling 
of a generic IR observable O for |0| > 1 (infinite slope at N") and |0| < 1 (vanishing slope at 

Nr)- 



Qualitatively, the scaling of physical observables near the conformal window shows a 
characteristic difference for |0| being larger or smaller than one. If |0| > 1 the value of a 
given observable as a function of the distance to the conformal window, O = 0{N^'^ — Nf), 
approaches the A^f axis with infinite slope. For |B| < 1 the slope vanishes on the quantum 
critical point A'^f = A'^^'^''. This characteristic shape dependence of the scaling is again a firm 
prediction of our scenario and may be observable in lattice simulations even away from the 
critical flavor number. 

Let us come back to the issue of choosing a specific fixing scale for comparing theories 
with different A^f to each other. As stressed above, e.g., the result for the shape of the 
phase boundary in Fig. |3| does depend on our choice of fixing the running coupling at the 
T-mass scale m^. This choice is not unique: in principle, the fixing scale can be chosen as 
a free function of A'^f. This would correspond to choosing an arbitrary function = kQ{N{) 
for the global scale occurring in the scaling relations (3.4) and ( [3.71 ). Of course, such a 
function induced by some ad-hoc scale fixing procedure could obscure our scaling relations. 
For instance, an extreme choice would be given by measuring all dimensionful scales in 
units of the critical temperature Tcr- In this case, the shape of the phase boundary would 
be a horizontal line at T/Tct = 1 terminating at A^f = Nf^ . Nevertheless, the scaling 
relations could still be verified, as they would translate into scaling relations for other 
external scales: e.g., the scale k at which the running coupling acquires a specific value 
(say a = 0.322) would diverge with A^'f — )• A^'^'^'^ according to /c ~ TcrlA'f — N^^\b for fixed 
Tcr (and < 0). This point of view can constitute a different way of verifying our scaling 
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relations on the lattice. In summary, these considerations demonstrate that the existence 
of scaling relations has a universal meaning, even though their concrete manifestation can 
depend on the details of the A'f-dependent scale fixing. In particular, the non-analytic 
structure governed by the exponent @ always remains. 



5. Scaling in low-energy models 

Let us study explicitly the scaling of two model observables, namely the quark condensate 
{qq) and the (constituent) quark mass Mg, by means of a simple ansatz for the effective 
potential U for the low-energy sector of QCD: 

Ui(t>) = -TT\n{d''+M^i^^)) + ^^Trln{d^+M^i(l)^))-NfN,TT\ni^+Mgi(l))) . (5.1) 

Here, (p represents a bi-fermionic scalar mean-field, the expectation value of which is related 
to the chiral condensate. Excitations on top of this condensate correspond to the Goldstone 
bosons, say the pions, and the sigma meson (radial mode). Mean-field expressions of this 
type generically arise in many low-energy QCD models such as NJL-type models or the 
quark- meson model. The Tr In terms simply correspond to the fluctuation contributions of 
the chiral mesonic and quark degrees of freedom. The masses and M„ of these mesons 
depend on dU/d(j)'^ and a linear combination of dU/dcp"^ and d'^U / dcfP'dcjP' , respectively, see 
e. g. Refs. |6|, 1^, g. The (constituent) quark mass Mg is given by the product of the 



Yukawa coupling h and (j). 

The expression Eq. (|5.lD for the effective potential is UV divergent and needs to be reg- 
ularized in some scheme (belonging to the definition of the model) at an effective regulator 
scale^ Agff . The IR observables can then be computed by solving Eq. (|5.1|) self-consistently 
for U{(j)) together with an initial renormalization condition provided at the scale Agfj, e.g., 
^ {'t') Ucft — l"^Aeff '?^^ + . . . . The bosonic mass parameter is naturally related to the effective 
scale, rwA^jj. ~ Aeg. 

For simplicity, let us consider the model in the limit of large number of colors, Nc ^ oo, 
for a fixed number of flavors Nf by way of example. In this limit only the quark loop 
survives and we are left with the following equation for the effective potential U: 

U{<P) = -iV/iVeTr In (^ + M,(<^)) . (5.2) 

Note that the Yukawa coupling is constant in this limit and can therefore be absorbed in a 
redefinition of the field (j), i.e., Mq = (j). The resulting equation for the effective potential can 
now be solved easily. The pion decay constant is given by the value of (j) which minimizes 
the potential U: 



dU 



= 0. (5.3) 

4>=4'o 



^Additionally, IR divergences can occur in the broken phase in the chiral limit, as the Goldstone bosons 
are massless. Such IR divergencies can be tamed, e.g., by the functional RG, yielding well-defined IR 
predictions. 
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We find 



(/-o ~ y^NjKAeS = ^/NjNlko\Ni - A^f r e , (5.4) 



where the last step holds near the conformal window, using the relation (3.5). Since the 
(constituent) quark mass is given by </>o, we conclude that Mq has a scaling behavior near 
NJ^ identical to that of the critical temperature. We would like to stress that the prefactor 
■sJNfNc is an outcome of our large- A'c analysis of the low-energy sector. In general, we 
expect that any observable O comes along with a complicated pre-factor function fo de- 
pending on the number of flavors Nf and N^- The determination of this function, e. g. for 
the constituent mass, may become complicated, depending on the truncations made in the 
low-energy sector. However, we stress that the A'^^ dependence coming from the prefactor 
is subleading compared to the scaling with [A'^f — A'^^^'l according to the IR critical exponent 
|0| of the running coupling. 

The scaling of the quark condensate can then be obtained by employing the following 
relation: 



\{qq)\ = ^ct>o. (5.5) 



Thus we find 



m)\-kl\N,-Nr\-^ . (5.6) 



The precise proportionality factor will again depend on the number of flavors and colors 
which, however, cannot modify the scaling behavior with respect to the distance to the 
conformal window. 

In ||4^, the quark mass, the chiral condensate and the pion decay constant have been 
computed within truncated Dyson-Schwinger equation for many-flavor QCD. Signatures 
of the quantum critical point have been identifled and the critical exponents have been 
extracted from a fit to the numerical data available away from the critical point for A'^f < 
Nf^ . In the light of our scaling relation, the results of [^9[ unfortunately remain somewhat 
inconclusive as A'^^'^'' has been fitted for each IR observable separately (yielding different 
values). This uncertainty is likely to spoil the fit for the critical exponents. We expect 
that a more careful analysis in the vicinity of the quantum critical point can easily put our 
scaling relation to test. 

Away from the chiral limit, the current quark mass is expected to modify the scaling 
relations. A generalized Gell-Mann-Oakes-Renner relation based on the fixed-point scenario 
in many-flavor QCD has been advocated in pO|. 



6. Conclusions 

Many-flavor deformations of real QCD are a fascinating testing ground for aspects of the 
chiral structure of QCD-like theories. The existence of a conformal window for larger flavor 
numbers A^f (below the critical value where asymptotic freedom is lost) gives rise to an 
interesting quantum critical point on the A'^f axis. Whereas a lot of effort has recently gone 
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into the determination of the value of A^f , we have concentrated in this work on the physics 
in the quantum critical region. Similarly to itinerant fermion systems [51|, we expect that 
the quantum critical point influences the properties of the phase diagram of a number of 
observables as a function of Nf in the neighborhood of the critical point. Note that the 
quantum phase transition in A'f can also be viewed as a higher-dimensional analogue of the 
Berezinskii-Kosterlitz-Thouless transition [52|. Similar transitions even occur in quantum 
mechanics which can also be analyzed in an RG language 

We have shown in this work, that the scaling of generic IR observables in the chirally 
broken phase as a function of N^^ — N{ exhibits a remarkably universal behavior. Our 
arguments were based on only few assumptions: the existence of an IR Caswell-Banks-Zaks 
fixed point in the running of the gauge coupling in the conformal window (which holds by 
construction if the conformal window exists), and the existence of a critical value of the 
gauge coupling for triggering chiral symmetry breaking. Whereas the first assumption is a 
universal statement, the latter assumption needs to be fulfilled only in specific RG schemes 
and for certain definitions of the coupling. In other schemes and coupling definitions, this 
assumption may translate into an assumption, for instance, on a sufficiently strong critical 
behavior of the quark-gluon vertex. 

We find that generic IR observables in the xSB phase scale with the distance to the 
conformal window in a characteristic fashion which is governed by only one independent 
critical exponent. This critical exponent is directly related to the corresponding critical 
exponent of the gauge coupling at the Caswell-Banks-Zaks fixed point. This relates uni- 
versal quantities with each other: chiral IR observables and the IR critical exponent. It 
establishes a quantitative connection between the chiral structure (Tcr, (ijip), fn, etc.) and 
the IR gauge dynamics quantified by the critical exponent B. Most importantly, it is a 
parameter free prediction following essentially from scaling arguments. 

Recently, the interest in the conformal window of many-fiavor QCD has been revived 
by technicolor scenarios |54, ^ for the Higgs sector of the standard model which have 
attracted renewed attention ||5^ . As these scenarios are largely motivated by the hierarchy 
problem of the standard model, it is important for some so-called walking models that the 
technicolor sector remains in the vicinity of the quantum critical point over a wide range of 
scales. The scaling laws should therefore directly apply to such models, even though each 
model, of course, has a fixed A'^f. 

It should be stressed that varying A'^f as a control parameter for a quantum phase tran- 
sition corresponds to comparing different theories with each other. Such a comparison is 
not unique for non-conformal theories but requires a specific choice of a dimensionful scale 
which is used as one and the same ruler for different theories. In the present work, we have 
argued that choosing a mid-momentum scale such as the r mass has many advantages. For 
small iVf, this choice has lead us to a simple estimate for the iVf dependence of the critical 
temperature, which is in very good agreement with lattice simulations. This type of argu- 
ment has already been successfully applied to the PNJL/PQM model, where implementing 
the linear Nf scaling for small A^f has helped adjusting the physical parameters, leading to 
a significant improvement of the thermodynamics properties of the model |37]. 

Coming back to the scaling relation near the conformal window, we are aware of the 
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fact that, for instance, the relation (|3.4| ) is difficult to test by lattice gauge theory: neither 
the fixed-point scenario in the deep IR nor large flavor numbers in the chiral limit are 
easily accessible. However, given the conceptual simplicity of the fixed-point scenario in 
combination with xSB, further lattice studies in this direction are certainly worthwhile. 
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